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We recalculate the color-Coulomb potential to one-loop order, under the assumption that the effect of the 
Gribov horizon is to make i) the transverse gluon propagator less singular; and ii) the color-Coulomb potential 
more singular, than their perturbative behavior in the low-momentum limit. As a first guess, the effect of the 
Gribov horizon is mimicked by introducing a transverse momentum-dependent gluon mass term, leading to a 
propagator of the Gribov form, with the prescription that the mass parameter should be adjusted to the unique 
value where the infrared behavior of the Coulomb potential is enhanced. We find that this procedure leads to a 
Coulomb potential rising asymptotically as a linear term modified by a logarithm. 



I. INTRODUCTION 

One of the early ideas regarding the confinement problem 
was that the confining force might come from one-gluon ex- 
change [1—3]. The suggestion was that a dressed gluon prop- 
agator, possibly combined with dressed quark-gluon vertices 
and arranged in ladder diagrams, would lead to a linear po- 
tential. Of course, the notion that the confining force can be 
entirely explained by ladder diagrams built from one-gluon 
exchange must nowadays be considered a little naive. There 
are many (related) problems with such a proposal, in particu- 
lar (i) the existence of a long-range color dipole field around 
static sources; (ii) long-range van der Waals forces which 
would then have to exist among hadrons; (iii) group represen- 
tation dependence (Casimir scaling) rather than N-ality de- 
pendence of the asymptotic string tension; and (iv) the ab- 
sence of color-electric flux tubes, not to mention the absence 
of string-like properties of such flux tubes, which have been 
convincingly seen in numerical simulations (cf. ref. [4] and 
references therein). Nevertheless, if it were possible to re- 
liably calculate the long-range behavior of, say, the color- 
Coulomb potential, then this information might be useful as 
an input into more sophisticated pictures, such as the gluon- 
chain model [5], where the problems just mentioned can be 
alleviated. Furthermore, the simple fact is that the instanta- 
neous color-Coulomb potential, is linearly confining. There is 
ample numerical evidence of this behavior [6-8]. 1 

This article is an an attempt to derive the long-range color 
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Coulomb potential analytically in Coulomb gauge. There 
have been a great many efforts in this direction over the years; 
refs. [9-17] is a partial list, see also refs. [18-21] for similar 
efforts in covariant gauges. Here we will focus on a simple 
one-loop perturbative calculation, modified minimally by cer- 
tain features associated with the Gribov horizon. 



II. GLUON PROPAGATORS AND THE GRIBOV 
HORIZON 

The potential energy of two static quarks in color represen- 
tation r is given in terms of the logarithm of a Wilson loop 
around a rectangular RxT contour 

V r (R) = - lim ^logW r (R,T) , (1) 

where W r (R,T) is the vacuum expectation value of the Wil- 
son loop. Let the side of length T be oriented in the time 
direction. For the purposes of the present article, the color 
Coulomb potential is defined by counting only the one-gluon 
exchange contribution to logW, and this is 

Vc(R) = f / -0zg 2 ND u (k,h = O)(l-e ik -*) , (2) 

where Df A (k) = 8 ab D^{k) is the 44-component of the gluon 
propagator in Coulomb gauge, C r is the quadratic Casimir 
in representation r, and N is the number of colors. The R- 
independentpart of this expression is the self-energy contribu- 
tion, which we will return to. It was shown by Zwanziger [22] 



cal simulations we would know that the color Coulomb potential must be 
confining, albeit not necessarily linear. 
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that g 2 D 44 (k) is a renormalization group invariant, and there- 
fore does not depend, e.g. in the context of dimensional reg- 
ularization, on the arbitrary scale fl. Vc{R) is the static quark 
potential which would be obtained if we approximate the log- 
arithm of a timelike Wilson loop expectation value by the one 
dressed-gluon exchange term; cf. [23-25]. Other treatments 
focus exclusively on the instantaneous part of D44, obtained in 
the £4 —> °° limit, but in this article we will include also non- 
instantaneous contributions to the potential, and this leads to 
setting k\ = 0. 

Let us define the renormalization-group invariant 



V(k) = -g 2 ND 44 (k,k 4 =0) 



(3) 



This quantity was computed to one loop long ago [26, 27] 
and the answer (at large |jfc|) is 



V(k)=- 



1 



Applying the one-loop result 

g 2 (li)N=- 



48,T- 



log- 



V QCD 



we obtain 



V(k) 



k 3^ ^g. 



(4) 



(5) 



(6) 



QCD 



which is indeed independent, to this one-loop order, of the 
scale /i introduced in dimensional regularization. 

However, the perturbative expansion is based on an implicit 
assumption that, apart from the gauge-fixing condition, the in- 
tegration over gauge fields is unrestricted; there is no cutoff, 
for example, in the amplitude of gauge field configurations 
contributing to the functional integral. But we have known 
for many years that this assumption is wrong. In the lattice 
formulation, in particular, it is known that if all gauge copies 
are included, then the sum over the Faddeev-Popov determi- 
nants of each copy will vanish. This means that the expecta- 
tion value of any gauge-invariant observable would take on the 
nonsensical value 0/0, as was first pointed out by Neuberger 
[28]. In the continuum it is also believed, since the seminal 
work of Gribov [29], that the functional integral should be 
restricted to a single gauge copy per gauge orbit, as in the 
proposed restriction to the fundamental modular region advo- 
cated by Zwanziger [22]. It seems difficult to implement such 
a restriction in practice. At a minimum we can ask that the 
functional integral be limited to the Gribov region, in which 
the lowest eigenvalue of the Faddeev-Popov operator is posi- 
tive semi-definite, and in fact this is achieved automatically by 
the gauge-fixing algorithms employed in lattice Monte-Carlo 
simulations, which find local minima of 



X fr=l 



(7) 



It is the fact that the gauge-fixed configurations are local min- 
ima, rather than just stationary points, which ensures that all 



eigenvalues of the Fadeev-Popov operator are positive. One 
can even go a little further. Since the lattice Monte Carlo pro- 
cedure will never generate more than one configuration per 
gauge orbit in the course of a finite simulation, an additional 
restriction to one configuration per orbit is, in some sense, su- 
perfluous. 2 

The limitation to the Gribov region has two expected conse- 
quences. The first, which is true in both Landau and Coulomb 
gauge, is that the Gribov horizon will impose a cutoff on the 
magnitude of quantum fluctuations of the transverse gluon 
field. This is easy to check in special cases. For example, one 
can construct a (lattice-regularized) plane wave of some fixed 
amplitude, and compute the low-lying eigenvalues of the lat- 
tice Faddeev-Popov operator. As the amplitude is increased, 
the lowest non-trivial eigenvalue Xq decreases, and eventually 
becomes negative. Configurations with amplitudes such that 
Ao < are to be excluded from the functional integration. 

Gribov [29] suggested that the restriction to the Gribov re- 
gion would result (in Landau gauge) in a gluon propagator of 
the form 



Df v (k) = 8 ab 5 MV - 



k^k v 



1 



k 2 J k 2 + 



(8) 



fr- 



aud this propagator clearly vanishes at k 2 — s- 0. Zwanziger 
[31] derived this form by adding to the action a term which 
was intended to implement the restriction to the Gribov re- 
gion. Gracey [18] has calculated the resulting static quark 
potential to one loop, in Landau gauge, which results from 
the Zwanziger action. This potential turns out to be non- 
confining. 

Lattice simulations, however, have rather decisively shown 
[32-34] that the Landau gauge gluon propagator has a finite 
non-zero limit at k 2 — > 0, as is the case for a massive propaga- 
tor, i.e. 



D?*(*) = 5"(5 MV - 



ab 



kuk v 



1 



k 2 



(9) 



Of course this form cannot be exactly right either; the gluon 
propagator cannot have a physical pole and must somewhere 
violate positivity. Various more complicated forms for the 
gluon propagator, which agree with (9) at low momenta, have 
been put forward, e.g. [35], [36]. 3 

The corresponding situation in Coulomb gauge is not so 
clear, at present. For the transverse gluon propagator at equal 
times, the Gribov-Zwanziger proposal is that 



Dff(k) = 8 ab ( ft. 



kikj 



1 



(10) 



2jk 



2 Of course, if one is interested in a particular selection of gauge copies, 
such as the fundamental modular region, or the "B-gauge" [30], then it 
is necessary to transform away from the gauge copies generated by the 
standard algorithms. 

3 Recently, Zwanziger has suggested a reason why the original proposal in 
ref. [31] might have failed, c.f. ref. [37]. Dudal et al. [36] have proposed 
a modification of the original Zwanziger action, to bring the result more in 
line with the lattice results. 
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and numerical calculations by the Tubingen group [38] seem 
to support this proposal. However, other recent calculations 
by Nakagawa et al. [39] on time-asymmetric lattices, while 
supporting a vanishing gluon propagator at k 2 — > 0, suggest 
a slower approach to zero than the Gribov-Zwanziger form. 
Nakagawa et al. conclude that larger lattices will be needed 
to settle the precise power falloff as k —> 0. In the absence of 
decisive lattice data on this point, we will here investigate the 
consequences of the Gribov-Zwanziger form (10) and also, 
for the purpose of contrast, a simple massive transverse prop- 
agator 



Dff(k) = 5 ab ( 8ij 



k(kj 



1 



k 2 - 



(11) 



Either form is obtained by the naive replacement, in the inte- 
gration over transverse gauge fields, 



DA 



I 



where 



DA'f exp 



M 2 (k) 



d 4 k 1 



{2k) 



4 2 M2 (*M"' 



-(k)Ar(-k) 



m 4 /k Gribo v prop ag ator 
m 2 massive propagator 



(12) 



(13) 



and where the subscript G on the left functional integral refers 
to the restriction to the Gribov region, with A"' r the renormal- 
ized transverse gauge field. The replacement is closely related 
to Zwanziger's suggestion [31], formulated in Landau gauge, 
that the restriction to the Gribov region could be implemented 
by adding an additional term to the action, and this addition 
includes a mass term with M 2 (k) = m 4 /k 2 . 

The second expected effect of the restriction to the Gribov 
region is special to Coulomb gauge. Coulomb gauge is a phys- 
ical gauge, and it has a Hamiltonian containing a non-local 
operator 



1 



-VD 



(-V 2 )- 



1 



-VD 



(14) 



involving two factors of the inverse Faddeev-Popov operator, 
which is responsible for the Coulomb potential. Evaluated 
for a configuration directly on the Gribov horizon, where the 
lowest F-P eigenvalue is zero, this quantity is singular. As 
Zwanziger has pointed out [22], we may expect that most con- 
figurations in the Gribov region are quite close to the hori- 
zon, for essentially the same reason that most of the volume 
of a sphere, in a large number of dimensions, is concentrated 
in the near vicinity of the surface. But configurations close 
to the Gribov horizon ought to have an enhanced density of 
near-zero eigenvalues, as compared to the spectral density of 
—V 2 , and a numerical study of configurations generated by 
lattice Monte Carlo bears this out [40]. 4 Thus, another effect 



of restricting configurations to the Gribov region should be an 
enhancement of the color-Coulomb potential in the infrared, 
assuming (as in the free theory) that the infrared behavior is 
associated with the low-lying eigenmodes of the F-P operator. 

Thus we are led to explore the consequences of the follow- 
ing two assumptions: first, that the restriction to the Gribov 
region can be approximately implemented, as in (12), by the 
simple addition of a momentum-dependent mass term, and, 
second, that the value of the mass parameter must be such that 
the infrared behavior of the Coulomb potential is enhanced 
beyond the usual l/k 2 behavior. 5 The way in which this 
could happen is illustrated by the following over-simplified 
scenario: The mass term will regularize the infrared behav- 
ior of loop integrals, and one might hope (ignoring inte- 
grations over Feynman parameters and so on) that the main 
effect is something like the replacement of log(Jfc 2 /A 2 ) by 
log((fc 2 + m 2 )/A 2 ) in eq. (6). Then, just by tuning m = A, 
the color-Coulomb potential at low momentum becomes 



V{k). 



1 



A^ 
~l*| 4 ' 



(15) 



much as in the old Richardson proposal [1]. We will now see 
how close we can come to realizing this scenario. 



III. ONE-LOOP INTEGRALS IN FIRST-ORDER 
FORMALISM 

The Coulomb potential is directly related to the 44 compo- 
nent of the gluon propagator. If we denote by 5" & n^ v the one- 
particle irreducible contribution to the Coulomb gauge gluon 
propagator, and noting that 114, = for i ^ 4, then the 44 com- 
ponent can be expanded, as usual, in a geometric series 



D44(*) = i(l+Il44(*)i- 



1 



1 



k 2 l-g 2 NYl(k) ' 



(16) 



where g 2 NTl(k) = rLi4(A:) jk 2 . We would then like to calculate 
Yl(k,k4 = 0) to one loop, with the restriction to the Gribov re- 
gion approximated by adding a mass term to the gauge-fixed 
action. Even at the one loop level, the loop integrals are com- 
plicated and non-covariant, and some are difficult to evaluate 
by standard formulas. It turns out to be much simpler to carry 
out the calculation in the first-order formulation, which is of- 
ten used when dealing with Yang-Mills theory quantized in 
Coulomb gauge (see in particular [22, 42^4]). 



4 It is interesting that removal of center vortices removes this enhancement, 
and pushes a typical configuration away from horizon. 



5 The prescription here is similar to that in ref. [41], where a dimensionful 
parameter in the gluon propagator was adjusted to the precise point where 
negative Faddeev-Popov eigenvalues disappear. 
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The starting point for the first-order formalism is the Eu- 
clidean partition function for Yang-Mills theory fixed to 
Coulomb gauge 



Z(J) = J DA^ 8 [V- A] det[. 



x exp 



(17) 



where ./# = — V -D is the Faddeev-Popov operator, and the 
color indices on the gauge field and field strength tensor are 
not written out explicitly, but are left implicit. One then intro- 
duces an Ej field via the identity. 



exp 



dTxFA 



jY / D£, exp 



(iEiF 0i 



1 



(18) 



The E field is split into a transverse and longitudinal piece 
Ej = E\ r — dj§, and then one integrates out the A 4 field, 
which generates a delta-function enforcing the Gauss Law 
constraint. This is followed by integration over the field, 
which eliminates both the Faddeev-Popov determinant and the 
Gauss Law delta function. The details of how this goes can be 
found, e.g., in ref. [43], and the result is 



Z[J] = 

DA'[ I DE\ r exp 



j^W-^iEf+^-igJ^ 



dtd 3 xd 3 y(p c + gJ 4 ) x .,K[x,y,t,A tr ] (p c + gj 4 ) y<t 



(19) 



where 



p a c (x) = -gf bc A b ' tr (x)E^(x) 



(20) 



and B" = jEtjkFjk is constructed from the transverse A-field. 
The non-local kernel, providing the Coulombic part of the 
Coulomb-gauge Hamiltonian, is 



K[ X , yi tX]^[^-\^ 2 )^- l ] a x b y 



(21) 



Then 



8 ab D 44 (x-y) 



g 2 Z 8Jl(x)8J b A (y) 



/=o 



= (K ab (x,y,A tr (x 4 )))8(x 4 -y 4 ) 

d 3 Zl K" c (x,zuA tr (x 4 ))p c (z u x 4 ) 

x / d 3 z 2 K bd (y,Z2 1 A"'(y4))p d (z2,y4) 



(22) 



The contribution to one loop is obtained by expanding 
K ab (x,y 7 A"'(x 4 )) up to second order in the coupling. Since 
the product pp inside the integrals over z 1 . is already sec- 
ond order, we can set K to its zeroth-order value in the inte- 
grand. The result is 



8" b D 44 (x-y) 



fab 



U +3g 2 r d f !fb jd\dh2(- s 



x.y 

2 jracd j?bef 



(A^(zi)A{(z 2 )) (^) z 



d 3 z\d 3 Z2 



-V 2 



_V 2 / viVM/'^v-zy/uvira 1 _ V 2 

/ .V.J; \ / Zl,Z2 

{{A c i {zux 4 )A e j {z2 7 y 4 )) (> {E d {zx,x 4 )E f j {z2,y 4 )) 



(dj) 



-V 2 



8{x 4 -y 4 ) 



-{A'i^x,x 4 )E f ] (z2 7 y 4 ))o{A){z2,y 4 )Ef{z u x 4 ))^ 



(23) 



Z2,y 



In ordinary perturbation theory, the zeroth-order propaga- 
tors are determined by simply removing the restriction to the 
Gribov region in the integral over A'[ . Introducing polariza- 
tion vectors 



with the usual properties 

he x (k) = , £*-*{k)ef{k) = 8 xx ' , (25) 



T i} ^e x *(k)ej{k) 
x 

j DAf (jfc) = / ' DA"{k,X) , (27) 
one can easily derive the zeroth-order momentum-space prop- 



A? r (k)=Z,e?(k)A a (k,X) 



(24) 



so that 
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agators in first-order formalism 



with 



(At(k)A b j (k'))o = 8 ab T i j(k) ¥ 8 4 (k + k') 

(Ef(k)E b (k')) = 8 ab T ij (k)^8\k+k') 
(E?(k)A b j (k!)) = S a %(k)^S 4 (k + k l ) 



(28) 



Taking eq. (23) to momentum space and inserting the propa- 
gators above, one finds for II (k) 



d 4 p Tijip) 



d 4 p Tij^Tjjjp-k) 
(2tt) 4 p 2 (p-k) 2 



) 4 P 2 (p-k) 2 

[p 2 -p 4 {p4-k 4 y 



(29) 



as originally obtained in ref. [43], see also [44]. The integrals 
can be evaluated under dimensional regularization, and the 
standard result for the one-loop momentum space Coulomb 
potential is obtained. 

Now suppose that instead of simply removing the restric- 
tion to the Gribov horizon in the integration over A ,r , we try 
to mimic its effect by insertion of a mass term, as in eq. (12). 
The effect on the zeroth-order propagators is readily obtained: 



(Af(k)A b j (k')) = 8 aO T i j(k) V2 



abn 



Ef(k)E^e))o = 8 a0 T iJ (k 



abn 



1 


k 2 - 


VM 2 {k) 


k 2 - 


fM 2 (ife) 


k 2 - 


hM 2 (k) 




k 4 


k 2 - 


^M 2 (k) 



8 4 (k + k') , 
8 4 {k + k') , 



The (unregulated) expression for Yl(k), in the £4 = case we 
consider here, then becomes 



U(k)=J l -J 2 



(31) 



where 



Ji =3 



kjkj r d 4 p 



Tijip) 



k 2 J {2n) 4 (p 2 + M 2 (k))(p-k) 2 

3£^ f_fp Tjj(p) 

2 k 2 J (2?r) 3 (p 2 +M 2 (k)y/ 2 (p-k) 2 



, (32) 



and 



h = 



d 4 p Tjj(p) Tjjjp-k) 
(2tt) 4 p 2 +M 2 (k) {p-k) 2 +M 2 {k) 
x[p 2 +M 2 (p)~p 2 ] 
1 r d 3 P mp-<Op- k Tij{p)Ti } {p-K) 



2k 2 J (2k) 3 (Op + (0 p - k (o p - k 



, (33) 



(o P =Jp 2 +M 2 (p) 



(34) 



Our task is to evaluate suitably regularized versions of J\ 2 for 
the two choices of M 2 (k) shown in eq. (13). 

IV. DIMENSIONAL REGULARIZATION, MASSIVE 
PROPAGATOR 

As a first step, we will compute the Coulomb potential to 
one loop using the massive transverse gluon propagator shown 
in eq. (11). We do not believe this propagator is correct in 
Coulomb gauge even at low momenta. In contrast to Lan- 
dau gauge, existing lattice simulations indicate an equal-times 
propagator which falls to zero at k 2 = 0, as already mentioned. 
The massive propagator is mainly useful as an illustration of 
how the potential can be enhanced by appropriately tuning the 
mass parameter, and also serves as a contrast to the results ob- 
tained in the next section. Technically, the massive propagator 
is simpler than the Gribov propagator case, in that standard di- 
mensional regularization can be applied without any difficulty 
to the relevant loop integrals. 

We now apply dimensional regularization, taking into ac- 
count the fact that & = 3 — 2e. Then 



U{k) = -{h-h)- 



2-2e 1 
— j2 — \ha + hb -h) + -jh , 



(35) 



where 



2f 



d 2m 'p 



(27t) 2co ' {p 2 +m 2 yi 2 (p-k) 2 ' 



kiki 
k 2 



ha 



hb = M 
I 4 = H 
7 5 = /i 



2f 



2f 



r d 2o) ' n 
2e 1 a P 




PiPi 


J (2n) 2<0 ' p 2 (p 2 + m 2 y/ 2 {p-k) 


■ d 2co p 


2 

m 




[2n) 2a (p 2 + 


m 2 )((p- 


-k) 2 +m 2 ) ' 


■ d 2m p 


P 2 




{2K) 2a) {p 2 + 


m 2 )((p- 


-k) 2 +m 2 ) ' 


■ d 2m p 


P\ 




{2n) 2a [p 2 + 


m 2 )((p- 


-k) 2 + m 2 ) ' 


■ d 2a p 


p^+m 2 


9 

-Pa 


(2n) 2m {p 2 + 


m 2 )((p- 


-k) 2 +m 2 ) ' 


p 2 k 2 -(p-k) 2 




P 2 (p-k) 2 







(36) 



and Co' = I — £, co = 2 — e. Integrals I\ through I 4 are diver- 
gent, I5 turns out to be finite. Before carrying out the usual MS 
subtractions, it is important to note that one is only allowed to 
make the subtractions which are made at m 2 = 0. In particular, 
one cannot subtract terms proportional to m 2 , because there is 
no counterterm which would generate such a subtraction. 

The integrals can all be evaluated by the standard methods, 
and the results for the divergent integrals are 



6 



1/1 \ 1 f /jfc 2 x(l-x)+m 2 x\ 

1/1 \ 1 /- 1 , , _i/2 /* 2 x 2 (l-x 2 )+m 2 x I \ 
72 = 125? U " 7+l0g4?r ) " 165?/o ^ " * - X2)Xi ^ j 

k 2 f x~ l,2 x 2 
+ ^ 2 Jdx l dx 2 9(1 - Xl -*2) k2x2(l l x2) 2 +m2xi ■ 

73 "=(^{G- r+1 ° g47r )-/^ l0g (^)}' 

7 - = ^ (I - r+i ° g47r ) - (^/"^ 2iog G?) - 1 (I - r+iog4?r+ (l +m2 ) 

74 = 4 (4iF G - r+ 1 +log47r ) (I + m2 ) + \jk 2 l dx Alog (£) • 



In these expressions we have defined 



A = k 2 x{l-x)+m z 



(38) 



All x-integrations run from to 1, and 0(x) is the Heaviside 
theta function. 

At this point we should take note of a source of possible 
trouble. In the first place, some of the integrals have produced 
m 2 /e terms, which cannot be subtracted away. Even finite 
terms proportional to m 2 would be catastrophic to our pro- 
gram, because these would tend to make the color Coulomb 
potential less, rather than more, divergent in the infrared. 
Somewhat remarkably, when the above integrals are inserted 
into (35), we find that there is a complete cancellation of 
the dangerous terms proportional to m 2 , while the remaining 
terms proportional to 1 /e — y+ log47T can be subtracted in the 
usual way. The end result is that 



n(*,M) 

3 



I6n 2 
3 



J dxx- 1/2 \og 



k 2 x{\ 



H ^ f dx\dx 2 d{\ — x\ — x 2 )x 

32n 2 J y ' 



-1/2 



X log 



k x 2 (l — x 2 ) +m 2 x\ 
~ J 2 



3 2 f 

k I dx\dx 2 Q[\ — x\ — x 2 ) 



I6n 2 
1 



k x 2 (l — x 2 ) + m 2 x\ 



i dxx(l — 2x)log i 
8tt 2 J \jj, 2 



1 



and therefore 



V(k) = -g 2 (fi)ND 44 = - 



48 n 2 k 



1 



k 



-n(t,/x) 



r/s , (39) 



(40) 



Inserting the one-loop expression for g 2 (/i), one finds that the 
dimensional regularization scale ji cancels out exactly, leav- 
ing the result 



V(k) = 



1 



(41) 



k z U(kA m ) 

Now we consider the infrared limit, k 2 /m 2 <C 1, starting 
with the integral 1$. Although this integral looks superficially 
divergent, it is clear, after an integration over p 4 which gives 

. 1 f d 3 p p 2 k 2 -(p-k) 2 



1 

2j (2tt) 3 p 2 (p-k) 2 ^(p-k) 2 + m 2 



y/p 2 +OT 2 - y/(p-A:) 2 +) 

^p 2 + m 2 +^(p-k) 2 +i 



(42) 



that in fact the integral is finite. Although it is still compli- 
cated, it is not hard to show that the low-momentum limit, up 
to 0(k 2 /m 2 ), is rather simple: 



1 



1 



360ft 2 m 2 



(43) 



It is also simple to evaluate the low-momentum limit of the 
single integrations over x: 

f wo / k 2 x(\ -x)+m 2 x\ 



/^x(l-2*)lo g 



4 r 

3 m 2 

2 



m 

■21og-5--4 

MS 



k z x(l 



A 2 - 

MS / 

1 k 2 1 m 2 

60 tn 6 AL- 
MS 



(44) 
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If this were all there were, then it would be possible to choose 
m o<: Aj^j so as to cancel the constant terms, leaving only a 
term proportional to k 2 /m 2 . This would lead to an overall 
1/Jfc 4 dependence for the color Coulomb potential, and there- 
fore to a linear potential. However, the integral I 2 leads to the 
two expressions involving integration over two Feynman pa- 
rameters, and these turn out to spoil the desired result. The 
double integrals can be evaluated analytically at low k 2 , with 
the help of the Mellin-Barnes transform and converse map- 
ping theorem [45, 46]. The details are reserved for Appendix 
A. The result, up to 0(k 2 /m 2 ), is 

f, , an .-1/2. ( k 2 x 2 (l-x 2 )+m 2 xi \ 
/ dx\dx 2 d(\-x\-x 2 )x x log I -J J 

V MS J 

4, m 2 32 n 2 f k 2 V /2 8 k 2 
= 3 l °Z^-- + -\ — I (45) 



MS 



4 \ m 2 ] 15 m 2 ' 



and 



/ dx\dx 2 Q(\ — x\ — X2)^j 

J krx?(l 



-l/2 y 2 
1 x 2 



x 2 (l —X2) +m 2 x\ 



3n 2 I k 2 \ 32 k 2 



8 \ TO 2 



15 m 2 



(46) 



Note the appearance of terms proportional to |Jk|. Therefore, 
at low momenta, 



V(k) 



1 



— I02- 

48 1U & 



1_ 
16 



( k 2 

64 \m? 



1/2 



151 k 2 
1440 m 1 



(47) 



We have suggested that m should be set to the unique value 
which would enhance the infrared behavior of the Coulomb 
potential. This value is 



= e 21 / 22 A 



MS ' 



leading to the final result at low momentum: 



64 



V(k) = -— e 21 / 22 -lMS 
{) 3 lifep 



(48) 



(49) 



Since the term proportional to |&| 3 is dominant at low mo- 
menta, this results in an asymptotic potential rising logarith- 
mically with quark separation. 



V. CUTOFF REGULATOR, GRIBOV PROPAGATOR 

The result found in the previous section would be a little 
disappointing, if the transverse gluon propagator actually had 
the massive form with M 2 (k) = m 2 . Tuning the mass parame- 
ter to the unique value which enhances the Coulomb potential 
does take us to a potential which rises faster than 1 /r, but the 
rise is still only logarithmic at large color charge separation. 
We will now investigate what happens in the (possibly) more 



realistic case where the transverse gluon propagator takes on 
the Gribov form. 

We again have Tl(k) = J\{k) — J 2 (k), where J\^ 2 are given 
in eqs. (32-34), but this time with the choice M 2 (p) = m 4 / p 2 . 
It is awkward to evaluate J 2 , in particular, by dimensional reg- 
ularization; one would end up with a complicated multiple 
integral over very many Feynman parameters. Since we are 
only interested in the small-A: 2 behavior of these integrals, we 
have found it convenient to follow a different strategy, based 
on a simple momentum cutoff at \p\ = A. 

We are aware that a momentum-cutoff regulator is danger- 
ous in gauge theories, and is likely to violate Ward identities 
and introduce spurious divergences, but these problems will 
not arise in our present one-loop calculation. This does not 
mean that the momentum cutoff procedure is necessarily con- 
sistent at higher loops, but that property is not crucial to us. 
What we are really after is to use the momentum cutoff re- 
sult to figure out what the one-loop result for V (k) would be 
in the MS scheme, without actually evaluating the integrals 
via dimensional regularization. This strategy requires that the 
momentum cutoff and dimensional regularization results can 
be matched exactly at one loop, by an appropriate choice of 
coupling g 2 in the cutoff regularization. That matching will 
be postponed to the next section. 

From this point on, since we will mainly be carrying out 
integration in three dimensions, we will denote 



*=l*l , P=\P\ 



(50) 



Of course the first equality is true even if k denotes the modu- 
lus of the 4-momentum, since we only consider the case where 
k 4 = Q. 

Begin with J\ , which, with a momentum cutoff, can be writ- 
ten as 

2 47T z Jo J -1 y/p4 -4- m 4 p- + k z - 2pku 

(51) 



and make the split 
1 



1 m 2 — \/ p 4 + m 4 



,4 m L 



m 2 \J p 4 + m 



4 ' 



so that 



h{k) = -{J\a+J\b) 



(52) 



(53) 



where 



1 f A 1 — u 2 

Ju= 4^Jo dPp3 L dU m 2 (p 2 + k 2 -2pku)> 



1 



1 m 2_ 



Jib = — — r / dpp / du- 

Ak 2 Jq J-i m 2 ^p 4 + m 4 

1-u 2 



p 2 +k 2 —2pku ^ 

Integral can be evaluated analytically, with the result up to 
0(k 2 /m 2 ) (and discarding terms of 0(1 /A 2 )) 



J\A — 



A 2 k 2 (ISlog § -46 



6n 2 m 2 m 2 4507T 2 



(55) 
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For the integral J\g we first expand in powers of k the term Then we write J2 = J2A + Jib where 



1-w 2 



p 2 + k 2 — 2p ku 



l-u 1 



2k (u 1 



■u) k 2 (-Au A + 5u 2 -i) , 
— + — 1 L + 0(k 3 ) 



(56) 

and find, again up to 0(k 2 /m 2 ) and discarding terms of 
0(1/ A 2 ), 



Jib 



1 



6QK 2 m 2 



2k 2 log 



2m 2 



10A 2 + 10m 2 log^ 



2A 2 



(57) 



Adding together J\a and J\b, we then have 
1 



/1 



3007T 2 



k 2 ( k 2 
15^r log^-log2 



46 — 



log ^ + log 2 



(58) 



Both and are quadratically divergent, but this is only 
an artifact of splitting J\ into two pieces. These quadratic di- 
vergences cancel exactly in the sum, as they must, since the J\ 
integral is only logarithmically divergent in the cutoff A. 

We employ a similar strategy to evaluate J 2 (k) at low mo- 
menta. Defining R p = y p 4 + m 4 , the integrand in (33) is 

\p-k\ \p-k\Rp-pR p - k 



F\(p,k,u) 



Rp-k \p-k\Rp + pR p - k 



1 1 



1 — u 



k 2 2p 2 +k 2 -2pku 



(59) 



Let Fo(p,k,u) be the same expression with R p ,R p _k both re- 
placed by to 2 

,2 



F (p,k,u) 



\p-k\ \p-k\-p ( 1 1 l-u 



p-k\+p\k 2 2p 2 +k 2 -2pku 



(60) 



1 f A 2 f 1 
J2A(k) = dpp J ^duF {) (p,k,u) , 

1 /-a ri 

1 ' • J- 



JiB{k) = ^ J dpp J ^du(Fi(p,k,u)-F (p,k,u)) . 

(61) 

The first integral can be done analytically, and again keeping 
terms to 0(k 2 /m 2 ) and dropping 0(1 /A 2 ), 



Ju(k) 



A 2 fe 2 (1051og^-3451og2 + 83) 



167T 2 OT 2 



72007T 2 /7i 2 



(62) 



To evaluate J2B, we expand the integrand in a power series 
in k. Then the integration over p and u can be carried out, with 
the result 



Jib = 



log; 



127 k 2 



1 A z 

16K 2 OT 2 



48 k 2 



Combining Ju and J2B 



Ji 



■log 



log-' 



7 



2A 2 



727T 2 



53 log 2 



(63) 



m 2 \ 9607T 2 6 m 2 2407T 2 9607T 2 



1 



48 7T 



, 2A2 



12k 2 



(64) 



As with J\ , the quadratically divergent terms in J2A and 
necessarily cancel in the sum, since the J2 integral is only log- 
arithmically divergent in A. 

Substituting the results for J\ and J2 into (31), we now have, 
up to 0(k 2 ), 



U(k) 



m 2 V 9607T 



and therefore 



41 k 2 
— log — 



73 



4007T 2 



log2 \ 
1927T 2 / 



11 2A 2 
— log — 



48 % 



12k 2 



(65) 



1 

7n 



41 



m 2 V 9607T 2 



k 2 

log — 



73 



log 2 



4007T 2 192a: 2 



11 



48 7T 



2A 2 
■log — 



12n 2 



(66) 



VI. CONNECTING THE REGULATORS 

Now that we have computed V(k) to one loop with a mo- 
mentum cutoff, the task is to figure out what the result would 



have to be in the MS scheme, because we would like to express 
our result in terms of a physical scale such as Aj^. The key is 
to show that it is possible to choose g 2 = g 2 (A) in the cutoff 
expression, such that an exact matching to MS is possible. 
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Denote 



II(*)= / [dp]R(p,k,m) 



(67) 



where [dp] denotes the multiple integration measure, and of 
course the integral is logarithmically divergent. What we 
would like to calculate is 

where the integral is dimensionally regulated, and the usual 
MS subtractions are carried out. What we actually compute, 
however, is 



1 



? 2 (A)N 7a 



- / [dp]R(p,k,m) , 



(69) 



where the integral is regulated with a momentum cutoff, and 
the dependence of the coupling on the cutoff is not yet speci- 
fied. This expression can be rewritten slightly as 

1 



mN -j A [dp]R { p, k ,0) 



where the second integration is finite, and needs no regulator. 
Now suppose it is possible to choose g 2 (A)N such that, as 

A — > oo 



1 



? 2 (A)JV J A 



[dp]R(p,k,0) 



1 



[dp]R(p,k,0) . (71) 



Then 



1 



L [dp]R{p ^ m) 

= 7k*-L [dp]R{p ' k ' 0) 

dp]{R(p,k,m)-R(p,k,0)} 



2 !„ w - L[dp]R(p,k,0) 
SjfgWN J ms 



1 



\dp]{R(p,k,m)-R(p,k,0)} 
- _[dp]R{p,k,m) . 



'MS 



(jl)N J MS 



(72) 



The conclusion is that if we can find a g 2 (A)N which satis- 
fies the matching condition (71) at m 2 = 0, then the cutoff- 
regulated calculation will give us the desired result in the MS- 
scheme for any m 2 . 
As before, 



V(k) = - 



1 



1 



(73) 



Starting with dimensional regularization and taking m 2 = 0, 
we have for J\ 

Jl = 48^ ( 12 (^-r+log47T-log-^)+28-241og2 



(74) 



while for J2, defining n = 3 — 2e, 

,2e r d n p 



h = n(n-l) 



m n -a> 



P-k 



El 

k 2 J (2n)"+ l (o p - k {(o p + (o p - k ) 
, 1 P 2 k 2 -(p-k) 2 ' 



n—l p 2 (p — k) 2 



(75) 



The second term in the squared parenthesis containing the 
combination ^yjfp leads to a convergent integral which 
can be done directly at n = 3 with the result 



48 n 2 



(16-241og2) 



(76) 



For the first term containing the unity, it is better to go back to 
D = 4 — 2e dimensions using the identity 



-J[dp]{R{p,k,m)-R{p,k,0)}, (70) JjPA 



p 2 -pI 



(O p - (0 p - k 



{Pi + (o 2 ) U + co 2 ) °>P-k(«>P + %-k) 



One then obtains 



1 



48* y(- e -r+^)-^j [ 2 + ^ 

Adding the two contributions, we find 

1 / , 1 . . k 2 , 53 



(77) 



(78) 



h 



48 7T 2 



- -7+log4^-log TT 2) + y -241og2 



and, altogether 



1/1 k 2 31 

Jl-J2 = ^o ( 11 (--r+log47T-log^) + — 



48 k 2 



(79) 



(80) 



When we compute the potential V(k) with a renormalized 
coupling g:Lg.(/i) in the MS scheme, then the terms propor- 
tional to £ — 7+ log47T can be dropped. 

Next we turn to the cutoff regulator. In this m 2 = case 



Jl 



d 3 p 



\-u l 



2 J (27r) 3 p(p 2 +k 2 -2pku) 



48^( 121 °4 +8 



(81) 



(82) 



where A is the momentum cutoff in this integral. For J2, the 
answer in cutoff regularization is 



Ji = 



1 f d 3 p (O p - (Op- k 



k 2 J (27T) 3 C0 p - k (C0 p + C0p- k 

A? 

487T 2 V"° k 2 



1 - 



IpV-(p-k) 



21 



2 p 2 (p-k) 2 



l — flog^ + 14-221og2 



(83) 
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Consequently, the final result is 



Defining A^g by the equation 



^1-/2 = ^2 (lllog^--6 + 221og2 



Now we equate the Coulomb potentials computed with mo- 
mentum cutoff and dimensional regularization at m 2 = 0. This 
means equating the denominators of (73), which is just the 
matching condition (71): 



11 



1 



1 



^(AW 48FV lll0g ^~ 6 + 221 ° g2 
1 1 / !£_ 31 



Therefore g 2 (A)N defined by 



(85) 



S 2 MS^ N 487r ' 



one finds that 



1 



g 2 (A)N 48 a: 2 



11 log 



log 



M 

AL- 
MS 



A 2 - 

MS 



49 
T 



-22 log 2 



(87) 



(88) 



is the choice of g 2 (A) required to convert our result from cut- 
off regularization to the MS scheme. 

A useful check of this method for converting cutoff 
regularization to the MS scheme is to go back to the massive 
propagator case in section IV, recalculate J\ and J2 with 
the cutoff regulator, and insert those values plus (88) into 
(73). When this is done, we find that our result agrees 
precisely with the result already obtained using dimensional 
regularization and MS subtraction, shown in eq. (47). 



g 2 m ^)N 48^ 



u 2 49 

Hl0gj=j + y-22l0g2 



_ (86) 
in terms of the MS running coupling 8jgc(l*)> is the cou- 
pling to be used to convert the momentum cutoff result to MS. 
I 



VII. THE COULOMB POTENTIAL, FINAL RESULT 



Inserting (88) into (66), we see that A cancels out, and the 



potential, in terms of the physical scale Ajjg, is 



v(*) = - ? 



11 

487T 2 " 



2m 2 59 

V MS 



1 k 2 

48 7t 2 m 2 



41 m 2 
20 1Og ^ 



log2 



219 
~25 



-1 -1 



(89) 



As in the case of the massive transverse propagator, we now 
set m 2 to the unique value at which power behavior of the 
Coulomb potential is enhanced in the infrared. This leads us 
to 

« = ^A M se 59 / 66 «1.73A Ki s, (90) 
and therefore, for k 2 <C m 2 

\%% 2 m 2 

y{k) = 7 : r , (91) 

k \2E l °gTT ~ — + ^5 ) 

where m is given in (90). We have finally ended up with a 
potential which behaves, in the infrared, as — 1 jk 4 modified 
by a logarithm. 

One often hears that a — 1 jk 4 potential in momentum space 
corresponds, upon Fourier transformation, to a linearly in- 
creasing potential in position space. Strictly speaking, this is 
untrue; the Fourier transform of — 1 /k 4 is actually minus infin- 
ity, due to the very singular behavior of 1 jk 4 as k — > 0. But this 
is precisely why it is important to include the quark-antiquark 



self-energies, as we have done in eq. (2). The Coulomb self- 
energies of quarks and antiquarks are also infinite, and this is 
not only the usual UV divergence which can be regulated with, 
e.g., a lattice cutoff. The Coulomb self-energies of quarks 
and antiquarks have, in addition, an infrared divergence, and 
a short-distance or high-momentum or lattice cutoff will not 
make this type of self-energy finite. In fact, this is already 
a reason why isolated quarks and antiquarks, or a non-singlet 
quark-antiquarkpair, are infinitely massive, and cannot appear 
as asymptotic states. But for a color singlet quark-antiquark 
pair, the infrared infinities of the self-energy and interaction 
terms precisely cancel, leaving only UV divergent contribu- 
tions to the self-energies, and a finite interaction term. This 
cancellation has been noted previously in ref. [47], in connec- 
tion with the instantaneous Coulomb interaction, where it was 
shown more generally that the cancellation of infinities is ex- 
act for any global color singlet combination of static quarks 
and antiquarks. 

The color Coulomb potential is 

Vc W = -%Jl0? vm - eik ' S) (92) 
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FIG. 1. Numerical calculation of — (k 2 V (k) ) _1 (solid line) compared 
to the infrared limit (upper dashed line) derived here, and the stan- 
dard one-loop perturbative result (lower dot-dash line). The x-axis is 
in units of in 2 . 



30 




2 4 6 8 10 

R 



FIG. 2. Potential V(R) vs. R, obtained from a Fourier transform to 
position space of the numerical solution for V(k). The result includes 
the self-energy, which is both ultraviolet and infrared divergent. The 
infrared divergence is cancelled, for a color singlet, by a correspond- 
ing term in the interaction, as explained in the text. To regulate the 
ultraviolet divergence we have made an arbitrary subtraction such 
that V(R) passes through zero at R = 1. V(R) is in units of m, R in 
units of 1/m. 

and, using the small-£ approximation (91) to V(k), the Fourier 
transform to position space gives us asymptotically 

, , R^oo C r ( 1207T m 2 \ 

^ = log(8.12, W */3)0 *- ^ 

This transform is carried out in Appendix B. However, the 
small-A: approximation is only valid at large distances, i.e. i? ^> 
1 jm, in which case the integral is sensitive mainly to the small 
k behavior of V(k). An expression for V(k) valid at all k will 
agree with (91) at small k, and the usual perturbative result (6) 
at large k. We do not have an analytical expression for V(k) 
valid at all k, but it is not hard to compute V (k) numerically, by 
evaluating J\ and Ji in eqs. (32-33) numerically. The result for 
(— k 2 V(k))~ l in cutoff regularization is shown in Fig. 1, and 
it interpolates nicely between our analytical result at small k, 
and the perturbative result for large k at m 2 = 0. 

The next step is to Fourier transform our result for V (k) to a 
potential in position space, i.e. eq. (92). Using the numerical 



result for V(k) at all momentum, of course the UV divergence 
of the self-energy will appear. On the lattice this UV diver- 
gence is regulated by the lattice spacing, and in the ordinary 
one-loop perturbative calculation of the Coulomb potential, 
the static quark self-energy is dropped altogether. In our case 
it is simplest to get rid of the UV self-energy divergence by 
making an arbitrary subtraction, such that the potential van- 
ishes aiR = 1 (in units of 1 jni); i.e. we compute V(R) — V(l). 
The result is shown in Fig. 2. 

It is interesting to compare our result with lattice data at 
large-A^. Of course one cannot directly compare string ten- 
sions, because of logarithmic modification of the linear term. 
The best one can do is to compare the slope of V (R) in Fig. 2, 
multiplied by the large-Af Casimir factor Cf/N = i, with the 
lattice result for the asymptotic string tension, extrapolated to 
large-A^. In order to make this comparison, we need the ex- 
pression for the lattice asymptotic string tension in units of 
Ajjj-. This has been derived in ref. [48] (see also [49]), which 
finds that Aj^/y/o = 0.503 (2) (40) at N -> <», where the un- 
certainties refer to statistical error, and an estimate of the sys- 
tematic error from all sources. Therefore the string tension 
at large-Af, derived from lattice Monte Carlo simulations, is 
a = 3.95A^. Since V C (R) = jV(R) doesn't really have an 
asymptotic string tension, the comparison with a depends on 
where we choose compute the slope of the Vq{R). At, e.g., 
R = 2mT x w 1.16A™ where a confining potential seems to 

MS 

have taken over from l/R behavior, we find 6 
«5.69A|- 

MS 

« 1.44(7. (94) 

VIII. CONCLUSIONS 

In this article we have explored the idea that, in Coulomb 
gauge, restriction to the Gribov region can be approximated 
by a momentum-dependent mass term in the action. Within 
the Gribov region, the bulk of configurations should lie near 
the horizon, and configurations near the horizon are expected 
to strengthen the long-range behavior of the color Coulomb 
potential. If the mass term should have this same effect, by 
suppressing (on average) configurations outside the Gribov 
horizon, then the mass parameter should be adjusted to the 
unique value at which the Coulomb potential is enhanced in 
the infrared. 

We have tested this idea at the one-loop level, by a per- 
turbative calculation of the non-instantaneous color Coulomb 
potential derived from g 2 Z)44 (fc, £4 = 0). For a momentum- 
independent mass term, the finding is that the infrared behav- 
ior is confining, but only marginally; the potential rises log- 
arithmically with quark-antiquark separation. However, for 



6 For finite N, a cml , at R = TmT 1 would be 5.69(1 - 1/W 2 )A^ 
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a momentum-dependent mass term leading to the propagator 
suggested by Gribov, the result is quite different: we find a 
confining potential rising as linear modified by a logarithm, 
and our potential is expressed in terms of the usual scale Ajjj-. 

This result, like most of its kind, must be interpreted with 
caution. In the first place, we have no idea how accurate our 
one-loop result may be. The best check would be to carry out 
the calculation further, to two loops, but this is a formidable 
task in Coulomb gauge. In the second place, we cannot be sure 
of the validity of the Gribov propagator in Coulomb gauge. At 
present the lattice Monte Carlo evidence is suggestive but not 
decisive on this point [38, 39], and we hope that our work will 
help to motivate further lattice investigations of this issue. Fi- 
nally, because of the logarithmic modification, the potential 
found here is certainly not an upper bound on the static quark 
potential. However, the upper bound derived by Zwanziger 
[9] only applies to the instantaneous Coulomb potential, rather 
than the full one-gluon exchange potential. It would be inter- 
esting to derive the instantaneous potential at one loop, along 
the lines we have followed here. Such a potential would be of 
particular interest for variational calculations of bound states 
and the gluon chain, and for this purpose the validity of the po- 
tential in an intermediate range of distances may be sufficient. 
We leave this case for future investigation. 
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is obtained by keeping all the singular terms in the Laurent 
series around each pole of M(f,s), in this case restricted to 
poles on the left-hand side of the fundamental strip. Then the 
Converse Mapping Theorem tells us that asymptotically, as 

x^O, 



p,k \ /" 



(A4) 



A proof of the converse mapping theorem is given in [45], and 
application to Feynman diagrams is found in [46]. 

The strategy is to put the integrals in eqs. (45), (46) in the 
form (A2), make a singular expansion of M(f,s), and apply 
the converse mapping theorm. For this purpose, we will need 
the Mellin-Barnes representation [50] 



1 



1 



(1+A) V Inih 



ds A 



_,i »r(v-j) 
r(v) 



(A5) 



The fundamental strip is in the region < Re(s) < v. 

Beginning with the integral (46), we apply the above iden- 
tity with v = 1 



7< 



/ = k / dxidx 



2 3/2 



1 f , (k l x 2 {\-x 2 ) 

x / ds ' 

2ni Jr 



m 2 x\ 



(A6) 



Interchanging orders of integration, the integrals over X\,x 2 
can be carried out exactly, with the result 



/ = 



2%i m 2 



ds 



% 2jnY(b-s) 
sm(ns) (2j-l)r(|-j) 



(A7) 



Appendix A: Application of the Mellin-Barnes transform 

The Mellin-Barnes transform M(f,s) of a function/, is de- 
fined by 



M(f,s)= / ,/.w' '/a; , 
Jo 

and the corresponding inverse transformation is 
1 



f(x) = —J r dsx- s M(f 1 s) . 



(Al) 



(A2) 



Now making a singular expansion, and applying the converse 
mapping theorem, we have 



1 k 2 f ( k 2 
I = t ds 



2%i m 2 Jr 



37Z 2 



£2\? 

2 

fll 



32 k 2 
15m 2 



3tt 2 /8 32 



64 



15s 105(j+1) 



(A8) 



The integral in (45) is handled in a similar way. First write 



The first integral is typically well-defined in a region 
(the "fundamental strip") of the complex s-plane, with 
Smin < Re(s) < Smax, and the contour T is a line parallel to the 
imaginary axis inside the fundamental strip. Now let 



P 



r p k 



k (s + p) k 



(A3) 



I = dxidx2&(l — x\ — x 2 )X} 



X lOE 



and use the identity 



m 2 x\ 

,4^ 



k 2 x 2 (l -x 2 ) 



m 2 x\ 



(A9) 



be a "singular expansion" (denoted "x") of M(f,s) on the 
left-hand side of the fundamental strip. A singular expansion 



1 f , n/s 

log(l+A ) = — / dsA -r-r — r 

2m Jr sm(ns) 



(A10) 
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where the fundamental strip is in the region —1 < Re(s) < 0. 
Then 



l ~ X2 , -1/2 I , m 2 

dx x x x I log-^-+logjC! 



/'= / dx 2 [ 
Jo Jo 

+ f dx 2 [ 
Jo Jo yjx\ 2ni Jr 



MS 



A ~ x i dx\ 1 



k 2 x 2 (l-x 2 )\ 5 n/s 



m 2 x\ 



sin Ks 
(All) 



Again interchanging orders of integration, carrying out the in- 
tegrations over x\,x 2 , and making a singular expansion, we 
find 



MS x 

njs 0fr(l-i) 



sin7Ti (1 +2s)T(j-s) 
4 m 2 32 If ( k 2 

— log — = — + ds\ 

5 * 2 9 Inih \m 2 



MS 



n 2 /4 8/15 



4 m 2 
3 l0g A^ 

MS 



i i + 1 J 

32 i?(#_y 



8 fc 2 



(A12) 



This completes the low-^: 2 evaluation of the integrals in (45), 
(46). 



Appendix B: Transform to position space 

In order to determine the asymptotic form of the one-loop 
potential as R — > °o, we will need to transform the momentum- 
space expression at small k, eq. (91), to position space. Ab 
sorbing a constant 
be written 



■ilog2- 



5^ into the logarithm, V(k) can 



V(k) 



960n 2 m 2 



1 



41 £4 lo „_*i 

g (8. I2n 

9607T 2 m 2 



41 



-Y(k) . 



(Bl) 



This is expected to yield a positive linear potential, modulo 
logarithms, plus an infinite constant which is removed by the 
self-energy term, as discussed earlier. 

Before proceeding, we should stress again that (Bl) is only 
valid at small k 2 <C m 2 . The excuse for taking the Fourier 
transform anyway is that the large-/? behavior we are inter- 
ested in is dominated by small k behavior, so the error at 
large k should only affect terms which are subleading in R. 
Note in particular that there is an unphysical Landau pole 
in (Bl) on the real axis, at a comparatively high momentum 
k = 8.12m = This pole is certainly not present in the 

result we have obtained numerically for V(k) at all momenta, 



which is displayed in Fig. 1. The Fourier transform of (Bl) 
will nonetheless require a prescription (e.g. principal value) 
for dealing with the unphysical pole, but the choice of pre- 
scription, as we will see, only introduces an ambiguity in sub- 
leading terms at large R. 

In the following we will switch to units m = 8.12m = 1, so 
that 



Y{k) 



k 4 \ogk 2 



(B2) 



The inverse log has a cut on the negative axis and a Lan- 
dau pole at k 2 = 1 . The discontinuity across the cut is easily 
evaluated, and the l/log£ 2 factor can be expressed through a 
Cauchy integral 



logk z Jo s + 



1 



k 2 k 2 -\ 



, with p (s) 



1 



log s + n 2 

(B3) 

In order to perform the Fourier transform, the IR singularity 
of the 1 /k 4 term is regularized by writing 



1 



1 



^ok 2 {k 2 - 



(B4) 



whose Fourier transform leads to an additional constant (infi- 
nite in the jj. — > limit), which is removed by the self-energy 
term. This removal amounts to subtracting V (0) from V(R). 

In the following we consider the dispersive (first term in the 
r.h.s. of (B3)) and Landau pole contributions separately 

f{R) = f{R) - y(o) 

= [V D (R)-V D (0)} + [V P (R) - Vp(0)} , (B5) 

where 



V D (R)-V D (0) 
d 3 k 



kR 



{2Kf 



l]V D {k) 



= lim 



d 5 k , , 



[e ikR -\ 



k 2 (k 2 +n 2 ) 



1 1 



log A: 2 k 2 -\ 



(B6) 



and 



Vp(R)-Vp(P) 
d 3 k 



I 



e' k R - l}V P (k) 



{2nf 



lim 

li^>0 



d 5 k 



JkR 



1] 



k 2 (k 2 + n 2 )k 2 -l 



(B7) 



14 



For the dispersive part one finds 



which after a few more manipulations can be written as 



Y(R) 



= lim 



d 5 k 



n4bj {2k f Jo 
d 3 k 

lim 



p 2 \k 2 k 2 + ii 2 Jk 2 



n^oJ (2n) 3 Jo 
11 11 



'rfsp (s) [«''*■* -1]^ 

jU 2 



1 1 



1 1 



s k 2 sk 2 + s ji 2 -sk 2 + n 2 ji 2 -sk 2 



lim / dsp(s) — 



2 \AkR 



+ - 



1 



R r ds 



n Jo s 



~ Rti -e~ R ^~ s ] 
l-e 



1+2- 



-(J?->0) 
R ^ 2 



SR 2 RyG 



Change variables s — >• z = sR , 

V D {R)-V D (0) 

R rdz 1 



in Jo z log 2 -jj + K 2 



1+2- 



(B8) 



Vb(J?) - Vb(0) 

R f m dz i 



?7T 7o z log 2 ^ + 7T 2 + 3 
R rdz 1 V? 



87TJO Z log 2 (z(!) ] 

R rdz 1 



^ Mog 2 (z(|) 2 )+^ 2 

R ndz 1 1 



JO Z l0E*U[ir I +JT 



R f l dz 



87T Stc/o Z log 2 fz(|) 2 U7T 2 ^+l 



(B9) 



Next, the term in the bracket is approximated by 



*xJi Mog 2 (z(|) 2 )+7r 2 V^+ 1 



(B13) 

The last integral is finite in the limit R — > °o, again leading to 
a term of the order of 0{R / \og 2 R). The remaining integral 
is dominated by z = and l/( v / z+ 1) can be expanded in 
powers of y/z leading to, in the limit R — > °o, 



dz 



1 



1 



nJo z iog 2 ( z (3) 2 )+^2 



R 



%n\og{R/3Y 



O 



log 2 R 



(B14) 



1+2- 



l-e"v^ 



2 Vi 



(BIO) 



which has the same limit in for both z — > and z — > °°, limits. 
Therefore in the integral 



■dz 



1 



z log 2 -^ + 7z; 2 



1+2- 



l_ e -v5 



v^+3 

(BID 

one can take the R — > °° limit since the resulting integral is 
convergent. This gives a contribution of 0(Rj log 2 R), which, 
as will be shown below, is subleading in the R °° limit, since 
the leading behavior is 0(R/\ogR). 

In the R — > °° limit, the leading behavior can therefore be 
obtained from 



V D (R)-V D (0)^-- 



R rdz 



nJo z log 2 ^- + 7T 2 y / z + 3 



, (B12) 



We now return to the pole term. For this we need to evaluate 



d 3 k 
p (2k)3 



JkR 



k 2 -l 



(B15) 



This integral is not well-defined, because there is a pole at 
k = 1 (or, in general units, k = 8.12m) on the positive real 
axis. The leading R dependence, however, does not depend on 
how the pole is circumvented. This is because, in the neigh- 
borhood of the pole, k is finite, while the leading-/? behavior 
is determined by the behavior of the integrand in the k —> 
limit. To illustrate this point, we consider a prescription "P" 
for how to skip the pole which excludes from the integration 
range the interval 1 — be < k < 1+ as, 



d k ik R 1 cos (R) + a sin (R) 



p (2nf 



k L -\ 



4nR 



(B16) 



where a = \og(a/b) parametrizes the ambiguity. For instance, 
for the principal-value prescription we have that a = b, and 
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thus a = 0. For arbitrary a, one finds In the limit of R — > °°, (B17) reduces to R/%71 which cancels 

the corresponding term in dispersive part, cf. (B13). So finally 
Vp(R) — Vp(0) me i eac ii n g behavior in the large-/? limit is given by 



,ikR 



87rlog(/?/3) 2 ^ \\og 2 R 



/• t/ 3 fc f ikR , l i_ 

~~ /l™ 7 (27r) 3 k 2 (k 2 + jU 2 ) k 2 — I or asymptotically, restoring constants and factors of m 

R / 2a „ 1 - (cos/? + a sin/?) \ R , . /1?f)7r m 2 \ 



(B17) 
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